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ABSTRACT 



In order to validate existing techniques for computer 
simulation of the welding process, an experiment is designed 
and carried out in which seam and girth welds are produced 
on a high strength steel cylinder. Fabrication of the HY-130 
cylinder, selection of welding parameters for the seam and 
girth welds, and the selection and arrangement of instrumenta- 
tion are described. The theoretical background of finite 
element computer simulation of heat flow in weldments as well 
as techniques for axisymmetric finite element analysis of 
thermal stress and metal movement during welding are dis- 
cussed. Finite element representation of weldment geometry 
for heat flow and stress analysis is presented. Metallurgical 
properties of the HY-130 cylinder seam weld are presented 
in the form of photo-micrographs and electron scan micro- 
graphs . 

Temperature, strain, and radial distortion are measured 
at various locations on the HY-130 cylinder and the measured 
values are compared in tabular and graphical form with those 
values predicted by use of the heat flow and stress analysis 
programs. Excellent agreement is achieved between measured 
and analytically predicted temperature distributions, but 
only qualitative agreement is shown for strains and radial 
distortion. 
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CHAPTER I 

I. A INTRODUCTION ; 

Welding has been used as the primary means of fabricating 
ocean structures for the past forty years. Design tools for 
the fabrication of welded structures have evolved primarily 
as the result of experimentation. These tools are usually 
expressed as empirical relationships among welding parameters, 
some of which include structure and weld joint geometry, 
physical properties of the structural materials and welding 
consummables , as well as methods and rates of energy trans- 
fer during the joining process. In most cases the use of 
these empirical relationships has been adequate to ensure the 
achievement of good quality weld joints in a production en- 
vironment. 

As design tools become more sophisticated, and the in- 
terrelationships among physical properties of welding and 
structural materials are better understood, design practices 
relating to allowable stresses in welded structures converge 
toward material limits. Development and refinement of tech- 
niques for detecting and classifying flaws in weldments is 
another design consideration which tends to increase confi- 
dence in welded structures subjected to more demanding re- 
quirements, allowing design engineers to increase allowable 
stresses or to reduce factors of safety or uncertainty. 
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An area of welding technology which requires further 
definition and refinement is the occurrence of residual 
stresses and structural distortion during and as the result 
of the welding process. Extreme non-uniform thermal cycles 
of weld metal and of structural metal in areas near the 
weldment produce unevenly distributed thermal stresses 
which cause localized plastic deformations, with resultant 
residual stresses and structural distortion. These residual 
stresses can significantly contribute to structural failure, 
particularly when the failure mechanism is dependent upon 
low applied stress, such as brittle fracture or failure re- 
sulting from stress corrosion cracking. Distortions caused 
by welding can also cause structural failure by reducing a 
structure's ability to withstand buckling loads. 

The ability to accurately predict the magnitude and 
type of residual stresses and distortion resulting from the 
welding of ocean structures would greatly enhance a design 
engineer's ability to estimate total stresses which materials 
must withstand during loading. More accurate estimation of 
total stresses would allow for more effective use of materials 
and the increased predictability of structural failure would 
ultimately result in the production of more reliable ocean 



structures . 
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I.B BACKGROUND ; 

Attempts to analytically predict or to simulate thermal 
stresses and metal movement which occur during the welding 
process have been the subject of considerable research effort 
during the past ten years. The complexity of the interrela- 
tionships among welding parameters which contribute to the 
formation of residual stresses in weldments does not allow 
the development of empirical relationships which are useful 
for more than one type of material, joint geometry, or welding 
process. Welding variables were too numerous and their in- 
terrelationships too complex for the development of a general 
technique for the analysis of transient thermal stresses, 
residual stresses, and distortion in weldments w Design 
engineers usually relied upon trial and error, or the use of 
design lanes derived from similar applications of proposed 
welding processes. 

Computer simulation of the welding process has offered 
reasonable opportunities to develop techniques which apply 
to wide ranges of materials, weld geometries, and processes. 

A description of the chronological and systematic development 
of techniques used for analysis of thermal stresses and metal 
movement during the welding process is contained in a paper 
presented by Masubuchi^ to an international conference on 
Residual Stresses in Welded Construction and their Effects. 
This conference was held in London 15-17 November 1977, 
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hosted by the Welding Institute, Abington Hall, Abington, 

Cambridge, England. Systematic development of computer 

techniques for welding simulation began with a one-dimen- 

2 

sional analysis by Tall in 1961. This research involved 

the simulation of thermal stresses occurring during the 

bead welding of a strip. Efforts at Massachusetts Institute 

of Technology by Masubuchi^ in 1968 expanded Tail's analysis 

into a one-dimensional FORTRAN program. This research was 

expanded into the development of a two-dimensional FORTRAN 

4 

program in 1970 by Iwaki . Iwaki attempted to apply two- 
dimensional finite analysis to bead on plate welding. The 
most recent computer program to be used at MIT was developed 

5 

by Muraki . Muraki's program uses finite element analysis 
techniques to predict thermal stresses and metal movement 
based upon temperature distributions within the weldment 
which are calculated by heat flow programs also developed 
at MIT. 

I.C SCOPE AND OBJECTIVES OF THIS STUDY ; 

The purpose of this research is to design and carry 
out an experimental procedure to generate data which can 
be used to validate the existing heat flow program and the 
program which uses finite element analysis to predict the 
occurrence of thermal stresses and distortion during welding. 
The production of seam and girth welds in a high strength 
steel cylinder under very carefully controlled conditions 
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while measuring temperature, strain, and distortion as a 
function of time and location with respect to arc travel 
would serve as a valuable aid in evaluating the accuracy 
and potential usefulness of the heat flow program and of 
Muraki's axisymmetric finite element program for analysis 
of thermal stress and metal movement during welding. This 
research is part of a continuing effort at MIT to validate 
computer analysis of welding. A summary of past experimental 
studies of themal stresses and metal movement during welding 
was presented by Masubuchi^ and is reproduced and updated 
here as Table I-A. 
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CHAPTER II 

MATERIAL CONSIDERATIONS 

II. A GENERAL ; 

The Office of Naval Research is sponsoring an ongoing 
program to investigate thermal stresses and distortion in 
high strength steel (HY-130) structures. HY-130 steel is a 
high quality, low carbon, quenched and tempered steel with 
a yield strength of 130-150 ksi. This material was produced 
by the United States Steel Company, primarily for use by the 
U . S . Navy . 

II. B COMPOSITION ; 

Acceptable ranges of chemical composition for HY-130 

steel are shown in Table II-A. These ranges, reported by 
20 

Shackleton, are composxtxonal Ixmxts of HY-130 steel as 
received for production. 

Mechanical properties of HY-130 steel in the "as received" 

21 

quenched and tempered condition, as reported by Schrodt, are: 

Yield Strength: 130-145 ksi 

Elongation in 2 inches: 15% 

Reduction of Area; 50% minimum transverse 

70% through thickness 

Impact requirements; 60 ft-lbf at 0*^ F 

A Rockwell hardness test conducted on the material to 

be used in this experimental study yielded a Rockwell "C" 

Scale Hardness of 32. 
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COMPOSITIONAL RANGES OF HY-130 STEEL (%) 

C 0.08 - 0.12 

Mn 0.60 - 0.90 

Si 0.20 - 0.35 

Ni 4.75 - 5.25 

Cr 0.40 - 0.70 

Mo 0.30 - 0.65 

V 0.05 - 0.10 

S 0.015 Max 

P 0.010 Max 

Ti 0.020 Max 

Cu 0.25 Max 



TABLE II-A 
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II. C TEMPERATURE DEPENDENCE OF MATERIAL PROPERTIES ; 

Subroutines within the axisymmetric finite element 
program developed by Muraki for analysis of thermal stresses 
and metal movement during welding require that material pro- 
perties be specified as a function of temperature. A search 
of the literature failed to produce detailed temperature 
relationships among the various material properties required 
by the Muraki program for HY-130 at elevated temperatures. 

Graphical presentation of these relationships reported by 
21 

Schrodt are reproduced in Figures II-l through II-IO. A 
tabular summary of this information is included as Tables II-B 
and II-C. 

The Muraki program which is to be verified by experi- 
mental data assiames a linear relationship between various 
material properties and temperature. Assumed linear relation- 
ships which served as input data to the Muraki program are 
indicated by a dashed line on the appropriate figures. 



STRESS (ksi) 
100 150 
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FIGURE II-l TYPICAL STRESS-STRAIN DIAGRAM FOR 
HY-130 AT ROOM TEMPERATURE 



YIELD STRESS (ksi) 
100 
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FIGURE II-2 ESTIMATED EFFECT OF TEMPERATURE ON 0.2% 
OFFSET YIELD STRESS FOR HY-130 



YOUNG'S MODULUS {ksi 
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FIGURE II- 3 ESTIMATED EFFECT OF TEMPERATURE ON 

young's modulus for HY-130 
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FIGURE II-4 TANGENT MODULUS CURVES FOR HY-130 
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FIGURE II-5 ESTIMATED EFFECT OF TEMPERATURE ON 
TANGENT MODULUS FOR HY^ISO 
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FIGURE II-6 ESTIMATED EFFECT OF TEMPERATURE ON 
POISSON'S RATIO FOR HY-130 
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FIGURE II-7 ESTIMATED COEFFICIENT OF THERMAL 
EXPANSION FOR HY-130 
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FIGURE II-8 EFFECT OF TEMPERATURE ON DENSITY OF HY-130 



THERMAL CONDUCTIVITY (BTU/hr ft F) 
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FIGURE II-9 ESTIMATED EFFECT OF TEMPERATURE ON 
THERMAL CONDUCTIVITY OF HY-130 



SPECIFIC HEAT (BTU/lbm F) 

.10 0.15 0.20 0.25 
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FIGURE II-IO ESTIMATED EFFECT OF TEMPERATURE ON 
SPECIFIC HEAT OF HY-130 
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SUMMARY OF MECHANICAL PROPERTIES 
FOR HY-130 STEEL 



TEMPERATURE 

(°F) 




0 


E 


H 


V 


68 


143 


30 


22.2 


.30 


200 


138 


29.6 


21.8 


.308 


400 


131.5 


28.8 


21.2 


.321 


600 


122 


27.7 


20.5 


.334 


800 


110 


26 


19.2 


.346 


1000 


88 


22.7 


17.0 


.359 


1200 


50 


14.0 


9.2 


.372 


1400 


25,5 


5.0 


4.9 


• .384 


1600 


11.2 


4.0 


3.1 


.397 


1800 


6.0 


3.0 


2.2 


.410 


2000 


4.0 


2.0 


1.5 


.422 


2200 


2.0 


1.2 


.8 


.435 


2400 


1.0 


.4 


.4 


.447 


2500 


0 


0 


0 


.454 



a = Yield Stress 2% Offset (ksi) 

-3 

E = Elastic Modulus (ksi x 10 ) 

_3 

H = Tangent Modulus (ksi x 10 ) 

V = Poisson's Ratio 



TABLE II-B 
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SUMMARY OF PHYSICAL PROPERTIES 
FOR HY-130 STEEL 



TEMPERATURE 

(OF) 




P 


'^T 


C 


K 


68 


0.284 


6.1 


19.4 


0.107 


200 


0.283 


6.5 


20.4 


0.11 


400 


0.282 


6.95 


21.5 


0.118 


600 


0.280 


7.3 


21.8 


0.126 


800 


0.279 


7.7 


21.2 


0.136 


1000 


0.278 


8.0 


19.5 


0.148 


1200 


0.276 


8.3 


17.1 


0.168 


1400 


0.275 


8.18 


15.0 


0.197 


1600 


0.275 


7.05 


13.4 


0.255 


1660 


0.275 


6.92 


13.1 


0.293 


1800 


0.275 


7.0 


14.0 


0.18 


2000 


0.275 


7.2 


15.0 


0.158 


2200 


0.275 


7.3 


16.0 


0.158 


2400 


0.275 


7.4 


16.7 


0.158 


2600 


0.275 


7.5 


17.0 


0.158 


2800 


0.275 


7.5 


7.0 


0.18 



3 

p = Density (Ibm/in ) 

Coefficient of Thermal Expansion 
(micros train/Op) 

C = Thermal Conductivity 

(BTU/hr ft °F) 

K = Specific Heat 

(BTU/lbm °F) 



TABLE II-C 
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CHAPTER III 

DESCRIPTION OF EXPERIMENTAL PROCEDURE 
III. A FABRICATION OF THE WELD SPECIMEN ; 

A circular cylinder of high strength steel was formed 
in order to test the MIT computer program for the analysis of 
heat flow in a weldment, with subsequent use of the tempera- 
ture distributions from the heat flow program to validate the 
Muraki axisymmetric finite element program for analysis of 
thermal stress and metal movement during welding. The HY-130 
circular cylinder was formed to dimensions suggested in the 
June 1977 proposal on Study of Residual Stresses and Distor- 
tion in Structural Weldments in High Strength Steels. This 
research proposal was made from Massachusetts Institute of 
Technology to the Office of Naval Research. 

Material for forming the cylinder was obtained from the 
Naval Ship Research and Development Center, Annapolis Labor- 
atory, Annapolis, Maryland. The material was received at 
Massachusetts Institute of Technology Welding Laboratory in 
the form of HY-130 plates 60" x 63" x 3/4". A request was 
made to Portsmouth Naval Shipyard in Portsmouth, New Hampshire, 
to roll the HY-130 plates into cylindrical shapes; Portsmouth 
Naval Shipyard being the only industrial activity in the 
vicinity of MIT with experience in forming HY-130. Detailed 
drawings and specific requirements were submitted to Portsmouth 
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Naval Shipyard and an agreement was reached for the shipyard 
to form the plates into cylinders. However, due to the 
thickness (3/4") of the plate to be formed to a small outside 
diameter (18") , Portsmouth Naval Shipyard engineers suggested 
that it -would be necessary to form a cylinder by bending two 
plates to semi-circular cylindrical shapes and then welding 
the two halves together. In view of this requirement and in 
order to facilitate easier handling of the HY-130 plates, the 
two plates were flame cut to four plates 60" x 31h" x 3/4" 
at the MIT Welding Laboratory prior to transporting to 
Portsmouth, New Hampshire. 

Stages of cylinder fabrication are depicted in Figures 

III-l and III-2. In order to prevent heat resulting from 

the flame cutting operation from affecting thermal properties 

22 

of the steel close to the proposed weld, Tsai suggested that 
flats required for bending be removed by flame cutting no 
closer than %" to the proposed leading edges of the semi- 
circular cylindrical mating surfaces. The flats were removed 
in this manner and remaining material to be removed in order 
to form a 60° Vee groove with 1/8" root, as illustrated in 
Figure III-2, was removed by mechanical grinding. The semi- 
circular cylinders were then transported back to the MIT 
Welding Laboratory for experimental welding. 
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ST/LGES OP CILDOER FABRICATION 



63 " ^ 




2i‘/t y /6 31^2" y 




FIGURE III-l 



stage 3 



Stage U 



Stags 5 



- Plate halves rolled to semi- 
circular cylindrical shapes 
TTith flats approximately 1 3/U” 
rKnaining on leading edges. 



- Flats removed ty flame cutting 
no closer than 1/U" to mating 
edges of semi-cylinders. 



- 1/U” excess material removed 

from mating edges of semi- 
cylinders and 30° angle pro- 
duced by mechanical grinding. 



FIGURE III-2 
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III.B SELECTION OF WELDING CONSUMABLES ; 

The type of welding suggested in the research proposal 
for girth welding of cylindrical shells was multi-pass gas 
metal arc welding CGMAW) . This type of welding requires the 
use of shielding gas and filler wire, both of which can 
significantly affect weld quality. The most common type of 
shielding gas used in GMAW of low carbon steel is Argon with 
25% CO 2 . However, in the case of HY-130, experience has 
shown that the use of Argon with 25% CO 2 leads to inadequate 

20 

fracture toughness properties in the weld deposit. Shackleton 
recommends rich Argon with less than 5% Oxygen or 10% CO 2 to 
avoid deterioration in the fracture toughness properties of 

♦ 

the weld deposit caused by Oxygen pickup from the shielding 
gas. Shielding gas selected and purchased for this experi- 
ment was Argon with 1% O 2 . This gas was reported by 
23 

Radziminski as the gas most effective in shielding HY-130 
welds . 

Filler wire selected for the HY-130 cylindrical seam 
and girth welds was 1/16" diameter Linde 140S wire. This 
wire was developed by Union Carbide specifically for HY-130 
weldments and was unanimously endorsed in the literature. Dur- 
ing periods when the filler wire was not in use, it was stored 
in a drying oven containing desiccative material. This step 
was viewed as necessary in order to prevent moisture or other 
impurities, or oxides, from collecting on the wire surface. 
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III.C SELECTION OF WELDING PARAMETERS : 

1. The variations in weldment geometiry which could 

be considered for this experiment were somewhat restricted 

by the fact that cylinder size required that the weld be 

made only from the outside. Based upon several previous 

24 

welding experiments reported by Doty and the above men- 
tioned restriction, a single 60° Vee groove was selected as 
near optimiam for producing the seam and girth welds of the 
cylinder. Weld joint geometry for this experiment is 
depicted below in Figure III-3. 




\ included 
angle 




8%" inside diameter 



1/8" mating surface 



FIGURE III-3 WELD JOINT GEOMETRY FOR PRODUCING 
SEAM WELDS ON HY-130 CYLINDER 
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2. For the purposes of this study, heat input (Q) is 
defined in the conventional way as the product of arc voltage 
and welding current divided by arc travel speed, as expressed 
below in Equation 3-1 : 

HEAT INPUT = Q = -■ (Equation 3-1) 

There were some variations in welding literature regard- 
ing the upper level of heat input which would produce a good 

20 

quality weld in 3/4" HY-130 steel. Shackleton recommended 

a maximum heat input of 40,000 Joules/inch. This upper limit 

25 26 23 

was also recommended by Nunez and ConnerT Radziminski 

recommended a maximvim of 45,000 Joules/inch. Although there 
was some variation in the suggested upper limits of heat 
input, there seemed to be general agreement that 40,000 Joules/ 
inch was a reasonable value for the production of good quality 
welds and this value was selected for production of the seam 
and girth welds on the HY-130 cylinder. Maximum and minimum 
allowable values of pre-heat and interpass temperature were 
also discussed in the above quoted references; general recom- 
mendations were more varied, but based on a rough average of 
recommended values, a common pre-heat and interpass tem.pera- 
ture of 200 was selected. 

3. Filler wire feed rate is automatically controlled 
by the G21AW welding machine . A wire feed rate of approxi- 
mately 87 inches/minute of the 1/16" Linde 140S filler wire 
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was measured during the production of a seam weld. The flow 

3 

meter for the Argon/1% 0 ^ shielding gas was set at 40 ft /hr. 
A sxammary of welding parameters which remained unchanged dur- 
ing the course of this experiment is shown in Table III-A 
below. The only parameters which were varied between welding 
passes were the arc voltage, welding current, and in one 
instance the arc travel speed. Voltage, current, and arc 
travel speed were recorded for each pass and are included in 
the presentation of results . 



WELDING PARAMETERS 


JOINT GEOMETRY 




60° Single V 
Groove 


MAXIMUM HEAT INPUT . . . 




40,000 Joules/in 


FILLER WIRE 




1/16" Diameter 
Linde 140S 


feed 


rate - 


87 in/min 


SHIELDING GAS 




Argon/1% ©2 


flow 


rate - 


40 ft^/hr 



TABLE III-A 
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III.D ARRANGEMENT OF APPARATUS : 

Welding equipment used for the entire experiment 
consisted of a Linde SVI-300 Power Supply with ST-2 Torch 
and associated arc travel and wire feed mechanisms. Seam 
welds were produced by placing the cylinder on an adjustable 
welding table and moving the welding torch along the groove 
of the seam weld. Prior to positioning the cylinder on the 
welding table, the two cylinder halves were tack welded at 
positions 6" from either end and at the center with 1" welds 
on the inside of each mating edge. Positions of the cylinder 
and welding torch for production of a seam weld are depicted 
in the photograph, Figure III-4, below. 




FIGURE III-4 POSITION OF CYLINDER AND WELDING TORCH FOR 
PRODUCTION OF SEAM WELD 
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III.E TESTING OF WELDING EQUIPMENT : 

In order to determine approximate temperature distri- 
butions during actual welding and to facilitate a thorough 
checkout of welding equipment prior to welding passes for 
data gathering, it was deemed advantageous to produce a seam 
weld on one side of the cylinder without benefit of instru- 
mentation to record temperature and strain. This allowed 
the wire feed mechanism, arc travel controls and indicators, 
gas flow regulators and meters, and the welding equipment 
electrical circuity and instrumentation to be tested and 
checked for accuracy under conditions expected during the 
actual test. This "test run" also provided an opportunity 
for operator training. A preliminary checkout of welding 
equipment had been conducted using "bead-on -plate" welds on 
scrap pieces of HY-130. 

Approximate temperature distributions to be expected 
during the instrumented test were determined by the applica- 
tion of temperature indicating liquid (TEMPILAQ) in narrow 
bands transverse to the welding direction. TEMPILAQ liquids 
corresponding to temperature levels of 400 °F, 600 °F, 800 °F, 
1000 °F, and 1100 ®F were applied to the chemically cleaned 
inner and outer surfaces of the cylinder. TEMPILAQ was re- 
applied between passes and distances of temperature penetra- 
tion from the weld centerline on both the inside and outside 
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of the cylinder were recorded for each welding pass . Welding 
current, arc voltage, arc travel speed, and maximum distances 
from the weld centerline for the various TEMPILAQ levels were 
recorded for each welding pass and are shown in Table III-B. 
This tabulated data was used as a basis for determining strain 
gage and thermocouple locations for the instrumented seam and 
girth welds to be used in validating heat flow and thermal 
stress computer programs . 
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AIDING PASS 


♦ TEMPERATURE PENETRATION DISTANCE 


1. I = 260A 
V = 23V 
T = 12 in/min 


OUKIDE INSIDE 

LEFT RIGHT LEFT RIGHT 

UOOF .86 - »9h .88 

600F - - .63 .60 

800F - - .UU .ho 

lOOOF - - .33 .35 

nOOF - - .25 .22 


2. I = 260A 

V = 23V 

V ■ 12 in/min 


UOOF .61 .38 .69 .75 

600F - - .38 .56 

800F - .22 .33 

lOOOF - - .06 .31 

llOOF - - - .19 


3. I = 260A 

V = 23V 

V * 12 in/min 


UOOF .80 - .88 .63 

600f - - .50 .38 

800F - - .13 - 

lOOOF - - - - 

HOOF - - - - 


U. la 26QA 

V a 23V 

V a 12 in/min 


UOOF - .92 .50 .83 

600F - .U9 .19 .50 

800F - - - - 

lOOOF - - - - 

HOOF - - _ - 


$. I a 3$0k 

V = 22V 

V a 10 in/min 


UOOF .92 1.55 1.06 1.63 

6oof - .80 .38 .75 

800f - .58 - .38 

lOOOF - .09 - - 

HOOF - .03 - - 


6. I - 33QA 

V a 21V 

V - 10 in/ min 


UOOF 1.U9 1.15 1.50 .88 

6 OOF .80 .75 .75 

800F .55 - .25 - 

lOOOF .U9 - - - 

HOOF .U9 - - - 


7. I = 35QA 
V a 22V 
T a 20 in/min 


UOOF .U6 .78 .38 .75 

600 F - .Uo - .13 

8 OOF - .25 - - 

lOOOF - .03 - - 

HOOF - - - - 


8. I a 370A 
V a 21V 
T a 20 in/min 


TEMPERATURE PENETRATION DISTAI4CES NOT RECORDED 


* Temperature penetration distances are expressed in inches, measured 
in the direction indicated from the weld centerline. 



TABLE III-B 
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III.F SELECTION AND ARRANGEMENT OF INSTRUMENTATION TO 

MEASURE STRESS AND DISTORTION DURING WELDING ; 

Stresses occurring during the thermal cycles associated 
with seam and girth welding of the HY-130 cylinder were de- 
termined by measuring strains as a function of time (i.e. - 
location of the strain measurement with respect to the time 
history of welding arc location) and as a function of distance 
from the weld centerline. Strains were measured by placing 
electrical-resistance wire gages on the outside surface of 
the cylinder. Bi-axial strain gages were placed with one 
axis aligned with the cylinder axis and welding direction 
(herein called the longitudinal direction) and the other axis 
aligned in a direction normal to the welding direction (herein 
referred to as the transverse direction. Characteristics of 
the strain gages used for this experiment are shown in 
Table III-C. 

In order to validate the MIT heat flow program, it was 
necessary to measure the heat distribution within the cylinder. 
Heat flow and temperature distribution within the HY-130 
cylinder was viewed as critical to validation of the Muraki 
program, inasmuch as strain, stress, and distortion calcula- 
tions are all based upon temperature distribution data which 
must be supplied as input to the program. Thermocouples used 
for measuring temperature distribution were Chromel/Alumel 
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STRAIN GAGE CHARACTERISTICS 


Type Used: 


FAET-12D-12S6 (EL) 


Manufacturer : 


BLH Electronics, Inc. 




Waltham, Massachusetts 02154 


Gage Factor: 


2.00 t .1% 


Resistance : 


120 i .2 Ohms 


Grid Length: 


0.125 inches 


Cement Used: 


BLH EPy-500 



TABLE III-C 



Strain indications were corrected for temperature 
by the following temperature error (^apparent strain) equation 

EAP = -91,37 + 2.17T - 2.33 x lO'^T^ + 5.31 x 10~^T 

-8 4 

- 2.40 X 10 T 
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wires which were peened into small (0.040") indentations 
drilled into the cylinder outer surface. Thermocouple 
circuit continuity was assured by the use of an inductive 
spot welder on the thermocouple connections . After circuit 
continuity was assured, the thermocouples were protected by 
the use of sealing cement. 

Since longitudinal and transverse strains were to be 
measured with gages which required temperature correction , 
it was necessary to place thermocouples in positions which 
would give temperature-time histories consistent with strain 
gage locations. Arrangement of strain gages and thermo- 
couples on the HY-130 cylinder outer surface is shown in 
Figure III-5. Arranging the strain gages and thermocouples 
in this manner allowed them to be used for measuring strain 
and temperature distributions during the production of the 
second seam weld (the first seam weld was made without in- 
strumentation) and then to be used again after the cylinder 
had been cut into sections and mating edges prepared for a 
girth weld. Arranging the strain gages and thermocouples 
at a 45° angle with respect to the seam weld and at equi- 
distant locations with respect to each other allowed the 
strain gages and thermocouples to give strain-temperature 
distributions which could be adjusted in time so as to be 



coincident . 
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FIGURE III-5 
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It should be pointed out that the practice of locating 
strain gages and thermocouples so as to require time adjust- 
ment of their readings assumes quasi-steady-state heat flow 
(i.e. - the heat flow/temperature distribution with respect 
to the welding arc is the same at a position 40” along the 
length of the cylinder as at a position 41" along the cylin- 
der) . This assumption is necessary in view of the fact that 
strain gages and thermocouples cannot be located at the same 
position. 

Arranging the strain gages and thermocouples at a 
45^ angle with respect to the seam weld was viewed as ad- 
vantageous since this scheme places the strain gages and , 
thermocouples which were subjected to the most severe thermal 
cycle during production of the seam weld in a position which 
is subjected to the least severe thermal cycle during the 
girth weld. Failures occurring during the seam weld are 
most likely to occur at locations which are least critical 
to measurements to be taken during the girth weld. 

After determining and marking the desired strain gage 
locations, strain gages were cemented to the cylinder outer 
surface. Physical size and weight of the cylinder prevented 
the usual technique of curing the strain gage bonding cement 
in laboratory ovens. In order to assure that the strain gage 
bonding would withstand expected temperatures, a 5 foot by 
5 foot thermal insulating blanket (graciously manufactured 
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for this purpose by Bethlehem Steel Corporation Shipyard in 
Boston, Massachusetts) was placed around the cylinder and 
electrical heating elements placed inside. The temperature 
of the cylinder was slowly raised to 250 and the strain 
gage bonding cement was cured at this temperature for a period 
of four hours, as directed by the strain gage manufacturer's 
instructions. The heating elements used for this purpose 
were later used to preheat the cylinder prior to welding. 

For production of the girth weld, it was necessary to 
cut the HY-130 cylinder into sections. Original plans called 
for cutting the 60 inch cylinder into two 30 inch sections 
which would have mating edges prepared for a girth weld by 
machining a 60° single Vee groove. Due to physical limitations 
of laboratory facilities it was necessary to reduce the length 
of the two sections to be joined by a girth weld to 20 inches 
each. This required that the HY-130 cylinder be cut into 
three sections, each 20 inches in length. This decision was 
made prior to production of the instrumented seam weld, and 
the strain gages were placed so as to permit measurements 
during the seam weld as well as during production of the girth 
weld. 

After completion of the seam weld, the 60 inch long 
HY-130 cylinder was transported to the Army Materials and 
Mechanics Research Center in Watertown, Massachusetts, where 
it was sawed into three sections, each 20 inches in length. 
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A large power hack saw was used to cut the cylinder into 
sections. The saw was used on slow speed and the blade and 
cylinder were continuously cooled by an oil and water cool- 
ing mixture during cutting. This procedure was viewed as 
necessary in order to prevent heat generated by the cutting 
action from affecting thermal properties of the mating edges 
to be joined by girth welding. 

After the two saw cuts were made, producing the desired 
three 20 inch long cylindrical sections, a h inch thick ring 
was sawed from the section which would not be used in the 
girth weld; again paying particular attention to preventing 
heat buildup in the material being sawed. A diagram of the 
resulting sections appears below as Figure III-6. 




FIGURE III-6 
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A cross-section of the previously produced seam weld 
was removed from the % inch thick ring and later ground, 
polished, and etched for analysis of the metallurgical proper 
ties of the seam weld. Details of this procedure and a dis- 
cussion of the results are contained in Chapter IV of this 
study . 

Mating edges of the two sections to be joined by girth 
welding were then prepared for welding by machining a 60*^ 
single Vee groove as shown below in Figure III-7. Heat 
buildup during machining was prevented by the use of a con- 
tinuously flowing oil and water cooling bath during machining 




FIGURE III-7 
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The two 20 inch cylindrical sections were transported 
back to the MIT welding laboratory where apparatus required 
for production of the girth weld was assembled. In order to 
control arc travel speed during production of the girth weld, 
it was necessary to design and manufacture a supporting struc- 
ture for the cylinder which would permit cylinder rotation at 
a controlled speed with respect to the fixed position of the 
GMAW torch. This was accomplished by producing a structure 
which was designed to support the two 20 inch cylinders from 
within. Details of this supporting structure are shown below 
in Figure III-8. 




FIGURE III-8 
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The rectangular base of the supporting structure was 
bolted to the base of a tilting turntable and the 20 inch 
long cylinder halves were fixed to the support by force 
exerted on the interior walls of the cylinder by the threaded 
dowels. The turntable was then tilted so as to place the 
now 40 inch long cylinder (with 60° Vee groove located in the 
center) in a horizontal position. The threaded dowels were 
adjusted in order to provide for constant clearance between 
the stationary GMAW torch and cylinder outer wall during 
cylinder rotation. Turntable controls were adjusted and 
calibrated in order to achieve desired arc travel speeds. 

A photograph of the apparatus in use during production of the 
girth weld appears below as Figure III-9. 




FIGURE III-9 
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In order to measure radial distortion during production 
of the girth weld, two Bourns Linear Motion Potentiometers 
were placed at the inner surface of the cylinder. Character- 
istics of the linear motion potentiometers are shown below 
in Table III-D. 



ACCEPTANCE TEST 


RECORD 


FOR LINEAR MOTION POTENTIOMETERS 


TESTS PERFORMED 


SPEC LIMITS 


MECHANICAL TRAVEL 


4 

, . 0 . 44 in max 


ACTUATION FORCE 


, . 4 oz min 


TOTAL ELEMENT RESISTANCE 


.. 4750 to 5250 Ohms 


INSULATION RESISTANCE 

NOISE OR CONTACT 


, . 50 Megohms min 

at 500 V DC 


RESISTANCE VARIATION 


. 100 ohms max 


END SETTINGS 


. 0 . 1 to 5.0% 

each end 


LINEARITY 


^ 1.0% VR max 



TABLE III-D 
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Displacement transducer circuitry and arrangement of 
linear motion potentiometers within the HY-130 cylinder are 
indicated in Figures III-IO and III-ll respectively. Prior 
to mounting the displacement transducers in position within 
the cylinder, a calibration test was conducted by deflecting 
the displacement transducers through measured distances and 
recording resultant voltages. A plot of this calibration 
test appears as Figure III-12. After completion of calibra- 
tion tests, the linear motion potentiometers were installed 
within the cylinder as indicated in Figure III-ll. Holes 
were drilled and threaded in the cylinder support mechanism 
and the linear motion potentiometer support rods were 
threaded and screwed into these holes at diametrically 
opposed locations within the cylinder. Potentiometers 
corresponding to Channels 1 and 2 of the recorder were 
located at positions 1/2 inch and 1 inch respectively in 
the cylinder axial direction from the weld centerline. 
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CHANNEL 1 DISPLACEMENT TRANSDUCER 
SENSING ELEMENT LOCATED h INCH 
FROM WELD CENTERLINE. INSULATED 
FROM CYLINDER INNER SURFACE 
BY CERAMIC BEAD. 



DISPLACEMENT TRANSDUCER SENSING E, 
MOUNTING RODS, SCREWED INTO 
CYLINDER SUPPORT. 



CHANNEL 2 DISPLACEMENT TRANSDUCER 
SENSING ELEMENT LOCATED 1 INCH FROM 
WELD CENTERLINE. ALSO INSULATED. 



CYLINDER SUPPORT M.ECHANISM 



FIGURE III-ll DISPLACEMENT TRANSDUCER SENSING ELEMENT 

LOCATIONS 



BOURNS LINIPOT 
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Strain and temperature were measured at the indicated 
locations during production of the instrumented seam weld; 
and strain, temperature, and distortion were measured during 
production of the girth weld. Equipment used to record the 
measured values included the following; 

Minneapolis-Honeywell H1508 24 Channel Visicorder 

Three Helland Universal Bridge Balancer Units, 

Model no. 82-6 (ser no's 1332,1307,1312) 

16 Channel Strain Gage Calibrator 
(one arm bridge adapter) 

Gould Brush 220 Two-Channel Recorder (distortion) 

General arrangement and operation of this equipment 
during production of the girth weld are indicated below in 
Figure III-13. 




FIGURE III-13 
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III.G PRODUCTION OF THE INSTRUI4ENTED SEAI4 AND GIRTH WELDS ; 

1. The Seam Weld : After the HY-130 cylinder had been 

suitably supported and aligned; instrumentation and other 
apparatus arranged, checked and calibrated, the instrumented 
seam weld was produced by making seven passes. Figure III-14 
below shows approximate dimensions and location of each pass. 
Values of weld parameters for each pass are tabulated in 
Table III-E. 




FIGURE III-14 
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2. The Girth Weld ; Production of the girth weld 
required careful positioning of the rotating HY-130 cylinder 
with respect to the GMA welding torch. The girth weld required 
six passes with approximate dimensions shown below in Figure 
III-15. Girth weld parameters are also included in Table 
III-E. 




FIGURE III-15 
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SEAM WELD PASSES 


GIRTH WELD PASSES 


WELD 
PASS NO. 


I 


V 


V 


I 


V 


V 


1. 


320 


26 


0.20 


270 


26 


0.20 


2. 


320 


25 


0.20 


290 


25.5 


0.20 


3. 


320 


25 


0.20 


270 


26 


0.20 


4. 


300 


26 


0.27 


285 


25.5 


0.20 


5. 


285 


25 


0.20 


280 


26 


0.20 


6. 


280 


25 


0.20 


290 


26 


0.20 


7. 


290 


25 


0.40 


- 


- 


- 



I = Welding Current in Amperes 

V = arc voltage in Volts 

V = arc travel speed in inches/sec 



TABLE III-E SEAM AND GIRTH WELD PASSES 
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CHAPTER IV 

METALLURGICAL CONSIDERATIONS 
IV. A PREPARATION OF SPECIMEN ; 

When the 60 inch cylinder was sawed into sections, as 
previously discussed in Chapter III, a cross-section of the 
instrumented weld was cut from the h inch ring. The weld 
specimen was cut to approximately 3/4 inch by 3/4 inch and 
marked for identification as shown in Figure IV-1. Mounting 
for polishing was in thermo-setting bakelite resin, pressed 
into a Ih inch diameter by 3/4 inch long cylindrical shape 
under hydraulic pressure of 4200 psi at 250 °F for 2 minutes 
using a Buehler press. 

The specimen was then rough polished by successive use 
of 60, 120, 180, 240, 320, and 400 grit 3M TRI-M-ITE Silicon 
Carbide paper. Fine polishing was accomplished by successive 
use of l.Oy, 0.3y, and 0.05y alumina polishing powder. Etching 
for photo-micrographs required about 15 seconds with 5% NITAL. 
The 500X photo-micrographs. Figures rv-2 through IV-5, were 
produced by the use of a LEITZ WETZLAR microscope with NPL 
50X/0.85 lens, using Polaroid type 52 POLAPAN film. Fine 
polishing was repeated and the specimen was etched using 1% 
NITAL prior to the production of 1200X and 6000X electron scan 
micrographs. These electron scan micrographs appear as Figures 
IV-6 through IV-1 3. Locations of the photomicrographs and the 
electron scan micrographs are indicated on Figure IV-1. 
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FIGURE IV- 1 DIAGRAM SHOWING LOCATIONS OF PHOTO-MICROGRAPHS 
AND ELECTRON SCAN MICROGRAPHS. 3/4" SAMPLE 
WAS CUT FROM THE SEAM WELD (SHOWING CROSS 
SECTION OF WELD - WITH NOTCH IN UPPER LEFT 
CORRESPONDING TO WELD CENTERLINE) . AFTER 
POLISHING AND ETCHING, 50 OX PHOTO-MICROGRAPHS 
AND 1200X/6000X ELECTRON SCAN MICROGRAPHS WERE 
PRODUCED AT THE LOCATIONS INDICATED. 
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FIG IV-2 PHOTOMICROGRAPH OF HY-130 

WELD METAL. 50 OX. 5% NITAL. 




FIG IV-3 PHOTOMICROGRAPH OF HY-130 

FUSION ZONE. 500X. 5% NITAL. 
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FIG IV- 4 PHOTOMICROGRAPH OF HY-130 
HAZ. 500X. 5% NITAL. 




FIG IV-5 PHOTOMICROGRAPH OF HY-130 

PARENT JIETAL. 50 OX. 5% NITAL. 




FIG IV- 6 ELECTRON SCAN MICROGRAPH OF 
HY-130 WELD METAL. 1200X. 
1% NITAL. 




FIG IV- 7 ELECTRON SCAN MICROGRAPH OF 
HY-130 WELD METAL. 6000X. 
1% NITAL. 




% 
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FIG IV- 8 ELECTRON SCAN MICROGRAPH OF 
HY-130 FUSION ZONE. 1200X. 
1% NITAL. 



FIG IV- 9 ELECTRON SCAN MICROGRAPH OF 
HY-130 FUSION ZONE. 6000X. 
1% NITAL. 




FIG IV-10 ELECTRON SCAN MICROGRAPH OF 

HY-130 HAZ. 1200X. 1% NITAL. 




FIG IV-11 ELECTRON SCAN MICROGRAPH OF 
HY-130 HAZ. 6000X. 1% NITAL 



r 



ft 
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FIG IV-12 ELECTRON SCAN MICROGRAPH OF 
HY-130 PARENT METAL. 120 OX. 
1% NITAL. 




FIG IV-13 ELECTRON SCAN MICROGRAPH OF 
HY-130 PARENT METAL. 6000X. 
1% NITAL. 
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CHAPTER V 

DESCRIPTION OF ANALYTICAL TOOLS AND PROCEDURES 
V.A INTRODUCTION ; 

Analytical procedures used in this study to predict 
temperature distribution, strain, stress, and distortion 
occurring during production of a weldment were developed at 
Massachusetts Institute of Technology. Both the heat flow 
program and the program used to analyze theinnal stress and 
metal movement use finite element techniques and are written 
in the FORTRAN-IV language. 

Discussion of these programs within this study will be 
limited to a brief synopsis of their theoretical background, 
presentation of assumptions and applicable equations upon 
which they are based, and a summary of welding parameters 
pertinent to production of the HY-130 weldment which were 
used as input to the programs. Detailed discussion of these 
programs, lists, card decks, and manuals describing their use 
are located within the Department of Ocean Engineering at 
Massachusetts Institute of Technology. 

V.B DESCRIPTION OF THE HEAT FLOW PROGRALl ; 

1. Theoretical Background . The finite element computer 

program used in this study to generate temperature distribution 

in the HY-130 cylinder girth weld is Program B-4 in a report 

27 

developed by Muraki and Masubuchi entitled "Computer Programs 



Useful for the Analysis of Heat Flow in Weldments" at MIT in 
1974. This program develops temperature distribution based 
on heat conduction in quasi-stationary state from a moving 
point heat source in a plate of finite thickness. The 
solution assumes; 

a. The size of the heat source is small compared 
to that of the body under consideration. 

b. Heat is supplied instantaneously. 

c. Thermophysical properties do not depend upon 
temperature . 

d. Heat is transferred in conduction only (other 
modes such as radiation and convection are not considered) . 

e. The material of the body is uniform and iso- 
tropic in the thermophysical properties. 

The solution is based upon the following equation which 
expresses the difference between the temperature at time t 
(in seconds after the arc is initiated) and the initial 
temperature : 



T-T = 
o 

Where ; 
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2 . Geometry of Weld Specimen and Description of 
Finite Element Representation of Temperature Distribution 
Program ; 

Division of the HY-130 cylinder girth weld into elements 
for use in the heat conduction program is illustrated in 
Figure V-1. Node coordinates (x, y, and z values) for each 
of the 114 nodes are listed in Table V-A. The weldment is 
represented by finite elements which are arranged so as to 
cause node locations for analytical determination of heat 
distribution to be coincident with locations of thermocouples 
which were actually placed on the cylinder outer surface in 
order to experimentally determine temperature distribution. 

Node coordinates listed in Table V-A are based upon an 
origin Cx=y=z=0) which is coincident with the starting point 
of the heat source (welding arc) , with subsequent times of 
interest corresponding to location of the heat source (in 
nodal coordinates) at the referenced time. This analysis 
assumes that the heat source moves on the plane z=0 along 
the x-axis. The x direction (direction of welding arc 
travel) would appear on Figure V-1 as normal to the y-z 
plane, positive direction into the paper. In this figure 
the y-axis is along the cylinder length and the z-axis is 
in the radial direction. 



2 (iNCHl^) FI NITE ELEMENT REPRESENTATION OF WELDMENT 

GEOMETRY FOR HEAT FLOW ANALYSIS 




FIGURE 



NODE COORDINATES (HEAT FLOW PROGRAM) 
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In order to determine temperature distribution at 
weldment nodes for the first pass of the HY-130 cylinder 
girth weld, the following parameters were used as input to 
the heat flow program: 

Velocity of welding arc (v) = 0.2000 in/sec 
Heat intensity (g) = 6084 Joules/^^^ 

Where; q = I ^ V x (Eq. 5-2) 

h 

I = welding current 
V = arc voltage 
h = plate thickness 
ri = welding arc efficiency 

d.2TC 

Density (p) = 0.2 82 lbm/j^^3 

Specific heat (k) = 124.4664 Joules/^j^^ Op 

Thermal conductivity (c) = 0.5250 Joule/^^^ Op 

It should be noted that this stage of the computer 
analysis does not provide for temperature dependence of 
material properties. Values used in this case to determine 
temperature distribution are those valid at 400 °F. It is 
also significant to note that for this analysis, an arc 
efficiency (r) ) of 0.65 was assumed. 
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No modifications were required to the analytical or 
mathematical techniques required for generating temperature 
distributions to be used by the Muraki strain, stress, and 
distortion calculating program. However, since large amounts 
of numerical data were generated by the heat flow program for 
input to the Muraki program, it was necessary to modify the 
heat program output. Since weldment geometry and finite 
element representation of the nodes and elements were the 
same (except for coordinates) for both programs, the heat flow 
program was modified so as to provide its output in the form 
of punched cards suitable for use as input to the Muraki 
program. This reduced the probability of introducing numeri- 
cal error betv/een the two programs as well as improving the 
ease of program interaction. 
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V.C DESCRIPTION OF THE PROGR?yyi TO CALCULATE STRAIN, STRESS , 

AND DISTORTION : 

1. Theoretical Background ; The program used in this 
study to analytically generate strain, stress, and distortion 
distributions within the HY-130 cylinder girth weld was de- 
veloped at MIT by Toyohiki Muraki for thermo elasto-plastic 
analyses of axisymmetric structural components. The Muraki 
program is written in FORTRAN IV language and uses finite 
element techniques in an incremental fashion to calculate 
strains, stresses, and distortion. 

The Muraki program has great flexibility in terms of 

allowing temperature dependency of material properties to be 

considered as well as considering the effect of melting and 

deposition of metal which occurs during welding. Weldment 

restraint and concentrated loads are also considered in the 

analytical solution. In the "Manual on Axisymmetric Finite 

Element Program for Analysis of Thermal Stress and Metal 

2 8 

Movement During Welding", Muraki presents a brief synopsis 
of the theoretical background upon which this analysis is 
based. This synopsis is reproduced here as follows: 

Principle of Virtual Work 

The principle of virtual work for the quasi-static 
problem is, in general, 

/// {6s}’^{a}dV - /// {6u}'^{P}dV - //_ {6u}'^{F}dS = 0 U) 

V V -1 
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where superscript T stands for the transpose of a vector or 

a matrix. {e}, {a}/ and {u} are vectors of the time rates of 

strain, stress and displacement, respectively. 

• » 

P and F denote the time rates of body force over the 
volume, V, and surface force on the boundary surface, 
respectively. 

The boundary condition on the specified displacement on 
the surface, 82 ^ is given as 

{u} = {u} on S 2 ( 2 ) 

It is noted that the problem here is formulated in the 
small displacement theory based on the rectangular cartesian 
coordinate system. 

Constitutive Equations 

One of the important features in the analysis of thermal 
stress and distortion during welding is the temperature de- 
pendency of material properties and the yield criterion. 
Therefore the constitutive equations are provided here in a 
general form. 

The rate of strain for the elasto-plastic state is given 

by; 

{e} = C3) 

where superscripts e and p refer to elastic and plastic com- 
ponents of the strain, respectively. 



83 



The rate of plastic strain is assumed as 

= A{g} (4) 

where A is a proportional constant and g is derivatives of 
the yield function, f, with respect to the stress, a. 

(i,j = 1,2,3); 



13 



rr - M 
^ 3C, . 



(5) 



Using Hooke's law and Equation (4), we have 

{e} = ^|De| ^ {a} + A{g} 
where ElDel is the matrix form of Hooke's law. 



( 6 ) 



Thermal strain and the temperature dependency of 
material properties are added to Equation (6) as follows; 



E|De| ^ {a}- I 2 - |De| ^{a} 
A + + A{g} 



(7) 



where e is the thermal strain, and 5. . is the Kronecker 

ID 

symbol . 



Equation (7) is rewritten for the rate of stress as. 



{a}= E|De| {c}+ - e® E|De|{6j^j} 



- AE|De| {g} 



( 8 ) 



* 

Note ; 



stands for the inverse of the matrix 
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If the von Mises' yield criterion is adopted, f in 
Equation (5) is given by 

f = 0 - C , T) (9) 

where C is related to the strain hardening and given as a 
function of plastic strain, and temperature, T. 

0 is referred to as the equivalent stress and given 
by 

0 ^ = 3 / 2 { 0 "}'^{ 0 "} ( 10 ) 

It is noted that { 0 ''} is the vector of deviator stresses 
which is different from { 0 }. 



When a state of stress in a portion of a body is plastic 
and then changes to a new state of stress which is still 
plastic, the following relation is obtained from Equation 
C9) : 



3f 

80 



0 



+ {• 



3f 



3e 



Ip) 






— T 
3T 



= 0 



( 11 ) 



Equation (11) with Equations (4) , (9) , and (10) leads 

to 

{0"}'^{0} = -2/30(- AH" + II T) (12) 

where 



H" 



3f 



o £ 

From Equations (8) and (12), A is obtained as follows: 



and 



1 
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A 



2/3 a H" + 

+ E 
E 



E{a"}'^| De| {g} 
{a} - £®E{a"}'^ 



(E{a''}'^ I De| {e} 
De|{6. j} 



+ 2/3 






( 13 ) 



Substituting A of Equation into (8) we have a general 
form of the stress-strain relationship as 

{a} = (E|De| - l/S^E|De|{g} E{a'' ] De | ) x 
({£} - £®{ 6 ^.}) 

+ E/E {|l| - 1/S E|De| {gHa"}'^) {ct} 

_ • ( 14 ) 

- 1/S^ E|De| {g} 2/3 ct T 

where 

S^ = 2/3 a H" + E{a"}'^|De| {g} 
and |l| is called the unit matrix. 

Equation (14) covers the stress-strain relationship in 
the elastic as well as plastic state. This can be easily 
confirmed by eliminating terms related to the yield function 
from Equation (14) . 

Yield Criterion 

As described in Equation (9) the yield function used 
here is that proposed by von Mises. 

The yield function is to separate the region of elastic 
deformation from that of plastic deformation. Since the 
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deformation is closely related to the loading process, the 
yield function is also called a loading surface which is 
defined in stress space. 

If a part of a material in question is in the region of 
plastic deformation, which means that the state of stress of 
the part is located on the loading surface, one of the 
following states will be realized at the next instance; 



where 



f = 0, f^ < 0 (unloading) 

f = 0, f^ = 0 (neutral) 

f = 0, f" > 0 (loading) 



f = 



3f 



3a 



■^ij + 



3T 



ID 



If a part of the material is in the region of elastic 
deformation, one of the following states will be realized 
at the next instance: 



f < 0 (elastic) 

f = 0 (plastic) 



(16) 



It is worth looking into Equation (15) more closely. 
Using Equations (10) and (12) with Equation (15) , the 
following relationship is obtained: 

f' = AH' (17) 

Equation (17) indicates that when non-strain hardening 
is assumed of a material. Equation (15) cannot be used to 
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judge whether or not a part of the material recovers from 
the plastic to the elastic state of stress. To avoid this 
difficulty, Equation (17) will be written in the form of: 



lim f' 



H "->0 



= A 



(18) 



H' 



(19) 



Therefore Equation (15) is written as, 

f = 0, A < 0 (unloading) 

f = 0, A = 0 (neutral) 

f = 0, A > 0 (loading) 

to cover non-strain hardening materials. 

Finite Element Approach 

Various types of axisymmetric finite elements can be 
derived based on Equation (1) . 

2 8 

In the computer program described in the Muraki 
Manual, a four node-isoparametric element is used. Follow- 
ing a standard procedure to derive an element stiffness 
matrix and element forces, the coordinates and the displace- 
ments for the elements are assumed by 



X = 



y = 



4 

Z h.x. 
i=l ^ ^ 






( 20 ) 



2 h.u. 
i=l ^ ^ 



u 



( 21 ) 




:'1 
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4 

V = 2 h . V . (21) 

i=l ^ ^ 



where h^, , h^, and are interpolation functions, 

and y. are nodal coordinates, and u. and v. are nodal 
1 11 

displacements . 

A key of the above derivation is use of the same in- 
terpolation functions for both the element coordinates and 
the element displacements. Details of the interpolation 
functions are 

h^ = %C1 - r) (1 - s) 

h 2 = %(1 + r) (1 - s) 

( 22 ) 

h 3 = %(1 + r) (1 + s) 

h^ = %(1 - r) (1 + s) 

where r and s are referred to as natural coordinates and vary 
from -1 to +1. 



The element strains in the axisymmetric condition are 




r in 'A 

3x 



3u 







Y 

3u ^ 3u 
3y 3x ^ 



( 23 ) 
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Because u and v in the strain are a function of x and 
y, the following relation is first required between the x-y 
and r-s systems 











9 




9x 


9y 


9r 


1 


9r 


9r 


- 


■ 


9x 


M- 




) 


9s 


9s 




(24) 



Derivatives of the displacements with regards to r and 



s are 



9r 



X 9h . 




u . 
1 



9s 



X 9h . 




u . 
1 



(25) 



9u 

9r 



X 9h . 




V. 

1 



u 

s 




V. 

1 



Using Equations (24) and (25) , the strains in Equation 
(23) are given in terms of nodal displacements, u^, u^, 
and u^ as follows: 



{£> 



B] {q} 



(26) 
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where 



r T-i-r 1 

{el ={e^. y^^I 



{q} = {u^, u^, u^, u^, v^, v^, v^, v^} 



and {B} is the matrix which gives a relationship between 
the strains and the nodal displacements. It is noted that 
{b} is a function of r and s. 



The stress-strain relationship of Equation (14) is 
rewritten in a simple form as : 

{a} = |d 1{£} + {a®} (27) 



where 




T 

xy 



} 



• 6 

{D} is the stress-strain relationship matrix and {cr } is the 
stress vector due to temperature change, and the temperature 
dependency of material properties and the yield function. 

The element stiffness matrix {k} is obtained using the 
principle of virtual work of Equation (1) and Castigliano ' s 
theorem as, 



{k} = ///{B}'^{DHB}y dedxdy 

2iT//{B}'^{DHB}y|j| drds 



( 28 ) 
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where 



3x 




3r 


3r 


3x 


3y 


3s 


3s 



Introducing a numerical integration, we have 

{k} = 2 ttS {B}T.{dHb} . . |j 1 . .Y. .a. . (29) 

•• Ij Ij Ijljlj 

-L , J 

where subscripts i and j indicate that ^ {B}, and [j] are 
evaluated at points i and j . ^ is weight factor at points 

i and j in the numerical integration. 

The program in the Muraki manual uses the Gauss-Legendre 
niimerical integration and in that scheme can be found in 

text books on numerical integration. 

The element force vector {f } corresponding to the second 
term in the right side of Equation (27) is obtained in a 
similar manner. 

{f®} = ///{B}"^{ct® Ivdedxdy 

T .ft 

= 2tt//{B} {a” }y I J I drds 

Again introducing the numerical integration we have 

{f®} = 2 ttS {B}T.{a®}Y. . 1 J| . .a. . (31) 

. ^ ID id' 'id ID 

-L t J 

In the current version of the computer program described 
in the Muraki manual, only the static load considered is the 
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is the concentrated load at nodal points. Therefore there 
is no other element force. 

The stiffness matrix and force vector for the whole 
body (i.e. the v/hole system) are obtained by assembling the 
element stiffness matrices and force vectors. 

The equilibrium equations for the whole system are 
obtained as 

E{k} {Q} + S{f®}. = Z{F}. (32) 

i i j J 

where i and j stand for the i-th element and the j-th node 
respectively . 

{Q} is the vector of the time rate of nodal displacement. 
{f} is the vector of the time rate of concentrated load. 

Equation (32) is further written in the incremental 
form as follows: 

E{k}.{AQ} + Z{Af®}. = E{AF}. (33) 

i ^ i ^ j ^ 



where 



AQ = QAt 

A jr® i ® A A. 

Af = f At 



AF = FAt 

At is the time increment. 



It is noted that the time increment is included in 



Equation (33) implicitly. 
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Once the nodal displacements for the system, {AQ}, are 
obtained, the element strains are calculated by Equation (26) 
and then the corresponding stresses by Equation (27) . 

When there is no further increment of concentrated and/ 
or thermal loads, the analysis is complete. 

When there is another increment of the loads , new 
equilibrium equations are calculated using the nodal dis- 
placements and stresses obtained in the preceding calculation. 
This procedure of incremental analysis continues until the 
final increment is solved. 
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2 . Geometry of Weld Specimen and Finite Element 
Representation of Weldment used by Muraki Program ; 

Finite element representation of the weldment is 
depicted in Figure V-2. A comparison of Figure V-2 with 
the finite element representation of the weldment. Figure V-1, 
used in the heat flow program will show that an axis trans- 
formation as well as a shift of the origin were required in 
order to use heat flow program output as input for the Muraki 
Program. 

For the Muraki program, the y-axis corresponds to the 
radial direction and the x-axis is along the cylinder length. 

The z-axis corresponds to the direction of welding arc travel, 
normal to the x-y plane as shown in Figure V-2, positive 
direction out of the paper. The origin (x=y=0) for the Muraki 
program analysis corresponds to a point at the center of the 
cylinder end. z=0 corresponds to the time of arc initiation. 

The Muraki program also required the specification of 
nodal coordinates (x,y values) . Element connectivity was input 
to the Muraki program by identification of the nodes, in counter- 
clockwise fashion beginning with the node closest to the heat 
source, associated with each element. Nodal coordinates input 
to the Muraki program are listed in Table V-B and element con- 
nectivity is shown in Table V-C. 

Thermophysical material properties required by the 
Muraki program were taken from Figures II-l through II-IO 
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in Chapter II of this study. Welding parameters required as 
input data for the first pass of the girth weld were the same 
as those used by the heat flow program. These parameters 
were discussed and values presented earlier in this chapter. 



Y(//^ch£S) finite element representation of weldment geometry 

FOR STRAIN, STRESS, AND DISTORTION ANALYSIS 
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CHAPTER VI 



Comparison of temperatures and values of strain, stress, 
and distortion during production of the HY-130 cylinder girth 
weld with those predicted by use of the heat flow and Muraki 
programs follows in the form of tabulated data and in graph- 
ical presentations. During production of the girth weld, arc 
travel began at the seam weld opposite the cylinder seam about 
which thermocouples and strain gages were placed. Times at 
which data comparisons are presented were selected as repre- 
senting a thermal cycle during which temperature, strains, 
stress, and distortion underwent greatest changes. 

Temperatures determined experimentally at the eight 
thermocouple locations for the times of interest are shown 
in Table VI-A. Distances from v/eld centerline are indicated 
beside the column headings designated Tl, T2, T3, T4, T4’, 

T3 • , T2 ' , and Tl ' . These temperatures are presented graph- 
ically in Figure VI-1. Temperature distribution predicted by 
use of the heat flow program is presented by a computer gener- 
ated graph. Figure VI-2. 

Qualitative comparison of the analytical and measured 
temperature distributions shows remarkable coincidence. The 
slightly higher overall values of the analytically generated 
values may be explained by the assumed welding arc efficiency 
of 0.65. 
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TABLE VI-A EXPERIMENTALLY DETERMINED TEMPERATURE DISTRIBUTION 
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FIGURE VI-1 GRAPH OF EXPERIMENTALLY DETERMINED TEMPERATURES 
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FIGURE VI-2 COMPUTER GENERATED GRAPH OF ANALYTICALLY PREDICTED TEMPERATURES 
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For the purpose of presenting comparative strain, stress, 
and distortion data, the following terros are defined: 



STR(*.*)X = EXPERIMENTALLY DETERMINED STRAIN 

(MICROSTRAIN) AT A POINT *.* INCHES 
FROM THE WELD CENTERLINE. 

STR(*.*)A = ANALYTICALLY PREDICTED STRAIN 

(MICROSTRAIN) AT A POINT *.* INCHES 
FROM THE WELD CENTERLINE. 

STSC*.*)X = EXPERIMENTALLY DETERMINED STPvESS (lbf/^^2) 
AT A POINT *.* INCHES FROM THE WELD 
CENTERLINE . 



STS(*.*)A = ANALYTICALLY PREDICTED STRESS (Ibf/. 2) 

in 

AT A POINT *.* INCHES FROM THE WELD 
CENTERLINE . 

DST(*.*)X = EXPERIMENTALLY MEASURED RADIAL DISTORTION 
IN INCHES AT A POINT *.* INCHES FROM THE 
WELD CENTERLINE. 

DST(*.*)A = ANALYTICALLY PREDICTED RADIAL DISTORTION 
IN INCHES AT A POINT *.* INCHES FROM THE 
WELD CENTERLINE. 



Positive values of strain and stress correspond to the 
tensile condition and negative values correspond to compressive 
strains and stresses. Positive values of radial distortion 
correspond to displacement in an outward direction from the 
cylinder center. Tables and graphs on the following pages 
use the terms as defined above to compare data. Strain and 
Stress directions correspond to directions shown in Figure V-2. 



f ^ 

i 
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Tables VI-B through VI-E present in tabular form 
comparisons between strains measured during production of 
the HY-130 cylinder girth weld and those predicted by use 
of the Muraki program. Tables VI-F through VI-I present 
the analytically predicted stresses. Experimentally deter- 
mined stresses were not tabulated since experimentally 
determined strains could not be transformed into stresses. 

There are three (ignoring creep) components of strain; elastic, 
plastic, and thermal. Elastic and thermal constitutive rela- 
tions are simple and the stresses resulting therefrom could 
be calculated if these strains could be separated into their 
components. However, the constitutive relationships among 
elastic, plastic, and thermal strains are complex, made further 
so by the presence of a strain hardening material and by the 
dependence of all material properties on temperature. For 
this reason, experimentally determined strains could not be 
transformed into stresses. Values of radial distortion measured 
during production of the girth weld and predicted by the Muraki 
program are shown in Table VI- J. Comparative graphs of measured 
and predicted data are shown in Figures VI-3 through VI-9. 
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LONGITUDINAL STRAIN 
(X-DIRECTION) 


TRANSVERSE STRAIN 
(Z-DIRECTION) 


TIME (sec) 


STR(1.0)X 


STR(1. 0) A 


STR(1.0)X 


STR(l.O) A 


125.0 


-44 


0 


200 


0 


127.0 


-130 


0 


235 


0 


128.0 


-153 


0 


277 


0 


129.0 


-177 


10 


318 


18 


129.5 


-190 


-2 


365 


32 


130.0 


-198 


-23 


370 


51 


130.5 


-198 


-127 


405 


132 


131.0 


-185 


-305 


420 


279 


132.0 


-185 


-420 


47 5 


391 


133.0 


-228 


-350 


607 


387 


134.0 


-390 


-271 


780 


401 


135.0 


-608 


-177 


922 


420 


136.0 


-823 


-68 


967 


446 


137.5 


-933 


127 


902 


481 


140.0 


-818 


474 


527 


531 


145.0 


-464 


1090 


-194 


600 


150.0 


-370 


1530 


-319 


641 


155.0 


-4 


1820 


-184 


668 


160.0 


0 


1930 


-129 


680 


170 . 0 


49 


2200 


-94 


714 


200.0 


176 


2250 


-49 


726 



TABLE VI-B EXPERIMENTAL AND ANALYTICAL STRAIN 1.0 IN FROM 
XVELD CENTERLINE 
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LONGITUDINAL STRAIN 
(X-DIRECTION) 


TRANSVERSE STRAIN 
(Z-DIRECTION 


TIME (sec) 


STR(1.5)X 


STR(1. 5) A 


STR(1. 5) X 


STR(1. 5)A 


125.0 


-245 


0 


320 


0 


127.0 


-310 


0 


365 


0 


128.0 


-325 


0 


395 


1 


129.0 


-345 


2 


410 


16 


129.5 


-345 


-7 0 


425 


27 


130.0 


-345 


-21 


450 


41 


130.5 


-340 


-105 


470 


92 


131.0 


-340 


-259 


485 


185 


132.0 


-330 


-358 


505 


260 


133.0 


-330 


-318 


570 


263 


134.0 


-320 


-297 


626 


277 


135.0 


-325 


-288 


681 


293 


136.0 


-330 


-283 


706 


311 


137.5 


-335 


-265 


750 


336 


140.0 


-340 


-205 


749 


373 


145.0 


-330 


11 


664 


433 


150.0 


-350 


276 


560 


477 


155.0 


-400 


527 


496 


513 


160.0 


1 -390 


1 651 


427 


530 


170.0 


-330 


1060 


257 


590 


200.0 


-291 


1490 


-246 


666 



TABLE VI-C EXPERII4ENTAL AND ANALYTICAL STRAIN 1.5 IN FROM 
WELD CENTERLINE 
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LONGITUDINAL STRAIN 
(X-DIRECTION) 


TRANSVERSE STRAIN 
(Z-DIRECTION) 


TIME (sec) 


STR(2.0)X 


STR(2 .0) A 


STR(2.0)X 


STR(2.0) A 


125.0 


-310 


0 


315 


0 


127.0 


-355 


0 


385 


0 


128.0 


-360 


0 


395 


1 


129.0 


-370 


-5 


410 


14 


129.5 


-385 


-16 


435 


22 


130.0 


-400 


-32 


440 


33 


130.5 


-395 


-115 


445 


68 


131.0 


-390 


-264 


455 


132 


132.0 


-380 


-365 


460 


186 


133.0 


-365 


-334 


475 


193 


134.0 


-355 


-321 


490 


207 


135.0 


-320 


-320 


498 


221 


136.0 


-310 


-325 


500 


236 


137.5 


-283 


-331 


502 


257 


140.0 


-225 


-337 


520 


288 


145.0 


-220 


-304 


555 


338 


150.0 


-384 


-216 


606 


378 


155.0 


-509 


-94 


666 


413 


160.0 


-636 


-16 


689 


431 


170.0 


-625 


304 


600 


496 


200.0 


-326 


836 


164 


587 



TABLE VI-D EXPERIMENTAL AND ANALYTICAL STRAIN 2.0 IN 
FROM D'TELD CENTERLINE 
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LONGITUDINAL STRAIN 
(X-DIRECTION) 


TRANSVERSE STRAIN 
(Z-DIRECTION) 


TIME (sec) 


STR(2. 5) X 


STR(2. 5) A 


STR(2.5)X 


STR(2.5) A 


125. 0 


-175 


0 


295 


0 


127.0 


-220 


0 


310 


0 


128.0 


-250 


0 


320 


1 


129.0 


-265 


-5 


325 


11 


129.5 


-265 


-8 


335 


16 


130.0 


-265 


-13 


355 


24 


130.5 


-265 


-30 


355 


47 


131.0 


-265 


-61 


365 


88 


132.0 


-265 


-86 


365 


124 


133.0 


-260 


-85 


365 


132 


134.0 


-260 


-89 


355 


143 


135.0 


-255 


-93 


355 


154 


136.0 


-245 


-98 


350 


166 


137.5 


-223 


-104 


340 


182 


140.0 


-125 


-114 


340 


205 


145.0 


-115 


-123 


355 


244 


! 150.0 


-249 


-118 


457 1 


276 


155.0 


-379 


-99 


552 


305 


160.0 


-501 


-83 


625 


320 


170.0 


-560 


1 


624 


377 


200.0 


-376 


203 


365 


463 



TABLE VI-E EXPERIMENTAL AND ANALYTICAL STRAIN 2.5 IN 
FROM WELD CENTERLINE 
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LONGITUDINAL STRESS 
(X-DIRECTION) 


TRANSVERSE STRESS 
(Z-DIRECTION) 


TIME ( sec) 


STS (1.0) X 


STS (1.0) A 


STS (l.O)X 


STS (1.0) A 


125.0 




0 




0 


127.0 




3 




22 


128.0 




143 




374 


12 9.0 




436 




593 


129.5 




207 




888 


130.0 




-295 




1236 


130.5 




-2658 




2669 


131.0 




-6614 




5369 


132.0 




-9464 




11431 


133.0 




-8670 




6405 


134.0 




-8263 




5400 


135.0 




-8235 




3860 


136.0 




-8420 




1908 


137.5 




-8619 




-1588 


140.0 




-8727 




-7665 


145.0 




-7337 




-16930 


150.0 




-4913 




-21859 


155.0 




-2501 




-24190 


160.0 




-1321 




-24732 


170.0 




2508 




-25156 


200.0 




6329 




-21790 



TABLE VI-F PREDICTED X AND Z DIRECTION STRESSES 1.0 IN 
FROM WELD CENTERLINE 
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LONGITUDINAL STRESS 
(X-DIRECTION) 


TRANSVERSE STRESS 
(Z-DIRECTION) 


TIME (sec) 


STS (1.5) X 


STS (1. 5) A 


STS(1.5)X 


STS (1.5) A 


125.0 




0 




0 


127.0 




0 




2 


128.0 




7 




29 


129.0 




222 




518 


129.5 




54 




741 


130.0 




-226 




1031 


130.5 




-2336 




1709 


131.0 




-6161 




2950 


132.0 




-8488 




4195 


133.0 




-7210 




4717 


134.0 




-6482 




5298 ■ 


135.0 




-6155 




5751 


136.0 




-6047 




6131 


137.5 




-5846 




6464 


140.0 




-5495 




6320 


145.0 




-4303 




4198 


150.0 




-2819 




1155 


155.0 




-1363 




-1724 


160.0 




-612 




-3085 


170.0 




2010 




-7334 


200.0 




5202 




-10709 



TABLE VI-G PREDICTED X AND Z DIRECTION STRESSES 1.5 IN 
FROM WELD CENTERLINE 



113 





LONGITUDINAL STRESS 
(X-DIRECTION) 


TRANSVERSE STRESS 
(Z-DIRECTION) 


TIME (sec) 


STS (2. 0)X 


STS(2.0)A 


STS (2 . 0) X 


STS(2.0)A 


125.0 




0 




0 


127.0 




0 




1 


128.0 




-6 




20 


129.0 




-13 




382 


129.5 




-293 




487 


130.0 




-716 




608 


130.5 




-3141 




612 


131.0 




-7513 




621 


132.0 




-10343 




974 


133.0 




-9197 




1626 


134.0 




-8632 




2215 


135.0 




-8436 




2661 


136.0 




-8476 




3045 


137.5 




-8508 




3563 


140.0 




-8571 




4237 


145.0 




-8168 




4976 


150.0 




-7258 




4859 


155.0 




-6154 




4135 


160.0 




-5498 




3548 


170.0 




-2939 




887 


200.0 




1115 




-3872 



TABLE VI-H PREDICTED X AND Z DIRECTION STRESSES 2.0 IN 
FROM WELD CENTERLINE 



114 





LONGITUDINAL STRESS 
(X-DIRECTION) 


TRANSVERSE STRESS 
(Z-DIRECTION) 


TIME (sec) 


STS (2.5) X 


STS(2.5)A 


STS (2.5) X 


STS(2.5)A 


125.0 




0 




0 


127.0 




0 




1 


128.0 




-2 




17 


129.0 




-23 




295 


129.5 




-50 




450 


130.0 




-89 




656 


130.5 




-269 




1301 


131.0 




-595 




2465 


132.0 




-827 




3487 


133.0 




-780 




3687 


134.0 




-772 




4001 


135.0 




-782 




4299 1 


136.0 




-807 




4618 


137.5 




-842 




5036 


140.0 




-889 




5658 


145.0 




-958 




6567 


150.0 




-1025 




7051 


155.0 




-1140 




7139 


160.0 




-1231 




7013 


170.0 




-1689 




5963 


i 200.0 
\ 




-2710 




2510 



TABLE VI-I PREDICTED X AND Z DIRECTION STRESSES 2.5 IN 
FROM WELD CENTERLINE 
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RADIAL DISTORTION 
(Y-DIRECTION) 


RADIAL DISTORTION 
(Y-DIRECTION) 


TIME (sec) 


DST(0.5)X 


DST (0. 5) A 


DST (1.0) X 


DST (1.0) A 


125.0 


0.0000 


0.0000 


0.0000 


0.0000 


127.0 


0.0000 


0.0000 


0.0000 


0.0000 


128.0 


0.0000 


0.0000 


0.0000 


0.0000 


129.0 


0.0000 


0.0000 


0.0000 


0.0000 


129.5 


0.0000 


0.0002 


0.0000 


0.0000 


130.0 


* 


0.0002 


* 


0.0005 


130.5 


* 


0.0009 


* 


0.0013 


131.0 


* 


0.0026 


* 


0.0028 


132.0 


0.001 


0.0035 


0.001 


0.0038 


133.0 


o.obi 


0.0029 


0.001 


0.0037 


134.0 


0.001 


0.0028 


0.001 


0.0037 


135.0 


0.001 


0.0028 


0.001 


0.0037 


136.0 


0.002 


0.0029 


0.002 


0.0038 


137.5 


0.002 


0.0032 


0.002 


0.0039 


140.0 


0.003 


0.0035 


0.003 


0.0041 


145.0 


0.003 


0.0039 


0.004 


0.0044 


150.0 


0.003 


0.0041 


0.004 


0.0045 


155.0 


0.003 


0.0043 


0.004 


0.0046 


160.0 


0.003 


0.0044 


0.004 


0.0047 


170.0 


0.002 


0.0046 


0.004 


0.0049 


200.0 


0.002 


0.0049 


0.003 


0.0052 



*DEFLECTION TOO SMALL TO QUANTIFY 

TABLE VI-J EXPERIMENTAL AND ANALYTICAL DISTORTION 0.5 IN AND 
1.0 IN FROM l-JELD CENTERLINE 
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FIGURE VI~3 COMPARATIVE PLOT OF EXPERIMENTAL AND ANALYTICAL STRAIN 1.0 IN FROM 
WELD CENTERLINE 
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FIGURE VI-4 COMPARATIVE PLOT OF EXPERIMENTAL AND ANALYTICAL STRAIN 
FROM WELD CENTERLINE 
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FIGURE VI -6 COMPARATIVE PLOT OF EXPERIMENTAL AND ANALYTICAL STRAIN 
FROM WELD CENTERLINE 
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FIGURE VI-7 PREDICTED X AND Z DIRECTION STRESSES 1.0 IN AND 1.5 IN FROM WELD 
CENTERLINE 



121 



00001 

( 



X 

1 


rv' 

1 


X 


rvj 

1 


1 

CE 


J 

cr 


1 

cr 


1 

cr 


O 


o 


LP 


LP 


PvJ 


C\J 


C\J 


C\J 


cn 


cn 


cn 


cn 


f— 


H- 


1 — 


f— 


cn 


cn 


cn 


cn 


□ 


0 


<] 


+ 




0059 0009 



005 - 



OOOh- 



005 ^- 



o 

OJ 



ooon - 



*NI/jgi] U(S*S)SiS QNb U(0'S]S1S 



FIGURE VI-8 PREDICTED X AND Z DIRECTION STRESSES 2.0 IN AND 2.5 IN FROM WELD 
CENTERLINE 
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FIGURE VI-9 COMPARATIVE PLOT OF EXPERIMENTAL AND ANALYTICAL RADIAL DISTORTION 
0.5 IN AND 1.0 IN FROM WELD CENTERLINE 
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CHAPTER VII 
CONCLUSIONS 

This study of residual stresses and distortion in high 
strength steel weldments has shown that existing computer 
programs require further refinement in order to be used as a 
day-to-day tool of the welding engineer. 

The occurrence of residual stresses and distortion during 
production of a weldment is a complex phenomena. Past efforts 
to simulate the welding process involved clever simplifying 
assumptions followed by attempts to carry out experiments which 
supported these assumptions. The new approach to welding 
simulation employs the enormous power of computer techniques 
to analyze v/hat actually takes place during the welding process. 
However, if a large computational process is used to model the 
complex welding process, the thermophysical material properties 
on which the analysis is based must be accurately known. 

The axisymmetric finite element program which was developed 
by Muraki for analysis of thermal stress and metal movement 
during welding is promising in that it showed good qualitative 
agreement with carefully generated experimental data. Heat 
flow and resulting temperature distributions within weldments 
can be accurately predicted by existing methods. Qualitative 
and quantitative comparison between analytical and measured 
tem.perature distributions showed rem.arkable coincidence. 
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CHAPTER VIII 
REC OMIEND AT IONS 

Inasmuch as analytical methods to predict stresses and 
distortion resulting from thermal cycles can only be as good 
as the thermophysical material properties on which they are 
based, it is recommended that studies be undertaken to 
accurately determine the temperature dependence of high 
strength steel properties at alevated temperatures. This 
would enhance further efforts to simulate the welding 
process . 

For further study of computer simulation of multi-pass 

r 

welding, it is recommended that a micro-processor be used 
to process welding parameter data. Reduction of strains 
requiring temperature compensation is a tedious task, 
especially if measurements are made at several locations 
and recorded on one photo tape. This additional source of 
error could be reduced or eliminated if parameters of 
interest could be calculated internally and displayed 
directly. 

High temperature strain gages could significantly con- 
tribute toward a more meaningful study. The 400 °F limitation 
of gages used in this study severely restricted their loca- 
tion. 
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